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Abstract. We study the Whitham equations for the defocusing complex modified KdV (mKdV) 
equation. These Whitham equations are quasilinear hyperbolic equations and they describe the 
averaged dynamics of the rapid oscillations which appear in the solution of the mKdV equation 
when the dispersive parameter is small. The oscillations are referred to as dispersive shocks. The 
Whitham equations for the mKdV equation are neither strictly hyperbolic nor genuinely nonlinear. 
We are interested in the solutions of the Whitham equations when the initial values are given by a 
step function. We also compare the results with those of the defocusing nonlinear Sclirodinger (NLS) 
equation. For the NLS equation, the Whitham equations are strictly hyperbolic and genuinely 
nonlinear. We show that the weak hypcrbolicity of the mKdV- Whitham equations is responsible for 
an additional structure in the dispersive shocks which has not been found in the NLS case. 
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1. Introduction. In [11, 12], Pierce and Tian studied the self-similar solutions of 
the Whitham equations which describe the zero dispersion limits of the KdV hierarchy. 

The main feature of the Whitham equations for the higher members of the hierarchy, 
of which the KdV equation is the first member, is that these Whitham equations 
are neither strictly hyperbolic nor genuinely nonlinear. This is in sharp contrast to 
the case of the KdV equation whose Whitham equations arc strictly hyperbolic and 
genuinely nonlinear . In this paper, we extend their studies to the case of the complex 
modified KdV equation, which is the second member of the defocusing nonlinear 
Schrodinger hierarchy The Whitham equations for the defocusing NLS equation are 
strictly hyperbolic and genuinely nonlinear, and they have been studied extensively 
(see for examples, [4, 6, 7, 10, 13]). However, for the mKdV equation, the Whitham 
equations arc neither strictly hyperbolic nor genuinely nonlinear. 

Let us begin with a brief description of the zero dispersion limit of the solution 
of the NLS equation 



Here Ao{x) and So{x) are real functions that are independent of e. Writing the solution 
tp{x, t; e) = A{x, t; e) exp ^v'^ S{x^t,e) ^ ^ using the notation p{x, t; e) = A'^{x, t; e), 

v{x,t;e) = dS{x,t,e)/dx, one obtains the conservation form of the defocusing NLS 
equation 




(1.1) 



with the initial data 




< 





dx"^ 



Inp 



(1.2) 
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The mass density p = and momentum density pv = ^^^(V'V'^ ~ i'*'4^x) have 
weak limits as e ^ [6] . These Umits satisfy a 2 x 2 system of hyperbolic equations 

dp d , , . 

f d (1-3) 

-ipv) + -iW + 2p^)=0, 

until its solution develops a shock. System (1.3) can be rewritten in the diagonal form 
for p ^ 0, 

where the Riemann invariants a and f3 are given by 

As a simple example, we consider the case with a = a =constant. System (1.4) 
reduces to a single equation 

f +2(a + 3/?)f =0. (1.6) 

The solution is given by the implicit form 

Pix,t) = f{x-2{a + 3l3)t) , 

where f{x) = /3{x,0) is the initial data for (3. One can easily see that if /3{x,0) 
decreases in some region, then I3{x,t) develops a shock in a finite time, i.e., dp/dx 
becomes singular. 

After the shock formation in the solution of (1.3) or (1.4), the weak limits are 
described by the NLS-Whitham equations, which can also be put in the Riemann 
invariant form [4, 6, 7, 10] 

-^ + Xg,i{ui,...,U2g+2)-g^=0 , i= 1,2,. ..,25 + 2 , (1.7) 

where Xg^i are expressed in terms of complete hyperelliptic integrals of genus g [8]. 
Here the number g is exactly the number of phases in the NLS oscillations with small 
dispersion. Accordingly, the zero phase g = corresponds to no oscillations, and 
single and higher phases g > 1 correspond to the NLS oscillations. System (1.4) is 
viewed as the zero phase Whitham equations. The solution of the Whitham equations 
(1.7) for g >1 then describes the averaged motion of the oscillations appearing in the 
solution of (1.1) (see e.g. [7]). 

Let us discuss the most important g = I case in more detail. We note that it 
is well known that the KdV oscillatory solution, in the single phase regime, can be 
approximately described by the KdV periodic solution when the dispersive parameter 
is small [1, 5, 16]. It is very possible to use the method of [1, 16] to show that the 
solution of the NLS equation (1.1) for small e can be approximately described, in the 
single phase regime, by the periodic solution of the NLS equation. The NLS periodic 
solution has the form 



p{x, t; e)= p3 + {p2 - ps) sn^{^ypl - pz e{x, t; e),s) 



(1.8) 
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with 9{x,t; e) = {x — Vit)/e and the velocity Vi = Vi{pi, p2, pa)- Here pi's are deter- 
mined by the equation obtained from (1.2) 



{P - Pl){P - P2){P - Ps) 



with pi > P2 > P3, and sn(z, s) is the Jacobi elhptic function with the modulus 
s = {p2 — P3)/{pi — Pz)- We can also write pi's as [3] 

Pl = + U2-U3- U4)^ 

P2 = - M2 + U3 - W4)^ (1-9) 



P3 = ^{Ui -U2-U3 + 

with Ui > U2 > Us > U4. The velocity Vi is then given by 

Vi = 2(ui + U2 + U3 + U4) . 

For constants Ui, U2, W3 and U4, formula (1.8) gives the well known elhptic solution 
of the NLS equation. To describe the solution p{x,t; e) of the NLS equation (1.2), the 
quantities Wi, U2, U3 and U4 are instead functions of x and t and they evolve according 
to the single phase Whitham equations (1.7) for g = 1. 

The weak limit of p{x,t;e) of NLS equation (1.1) as e — > can be expressed in 
terms of pi, p2, ps and p4 [6] 



E{s) 



p{x, t)= Pl- (pi - ps) , (1.10) 

where K{s) and E{s) are the complete elliptic integrals of the first and second kind, 
respectively. This weak limit can also be viewed as the average value of the periodic 
solution p{x,t;e) of (1.8) over its period L = 2eK{s)/^pi — p^. 

In order to see how a single phase Whitham solution appears, we consider the 
following step initial data for system (1.4) 



a{x,0) = a, /3{x,0) = { ' , (1.11) 




where a > b, a > c, b ^ c. The solution of (1.4) develops a shock if and only if 6 > c 
(cf. (1.6)). After the formation of a shock, the Whitham equations (1.7) with g = 1 
kick in. For instance, we consider the Whitham equations with the initial data [7] 



ui{x,0) = a, U2{x,0) = b, M3(a;,0) = < ' , U4{x,0) = c. (1-12) 




Now notice that the Whitham equations (1.7) for g = 1 with the initial data (1.12) 
can be reduced to a single equation u^t + Ai,3(a, b, U3, c)u3x = 0. The equation has a 
global self-similar solution, which is implicitly given hy x/t = A3(a, 6, U3, c). The x-t 
plane is then divided into three parts 

(1)|<7, (2) 7 < I < 2a + 46+ 2c , (3) | > 2a + 46 + 2c , 

where 7 = 2(a + 6 + 2c) - 8(a - c){b - c)/{a + b - 2c) (see (2.7) and (2.8) below for 
the derivation). The solution of system (1.4) occupies the first and third parts, i.e.. 
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(1) for x/t < 7, 



(3) for x/t> 2a + Ah + 2c, 



a{x, t) = a, (3{x, t) 



a{x, t) — a, (3{x, t) ^ c 



The Whitham solution of (1.7) with g — \ hves in the second part, i.e. 
(2) for 7 < x/t < 2a + 46 + 2c, 



ui{x, t) = a , U2{x, t) 



t 



Ai,3(a, 6, U3,c) , Ui{x,t) = c, 



where the solution can be obtained as a function of the self-similarity variable x/t, 
if 



9W3 



(a, 6,U3,c) ^ 0. 



Indeed, it has been shown that the Whitham equations (1.7) are genuinely nonlinear 
[6, 7], i.e.. 



du,. 



{ui,U2,U3,U4:) > 0, 1 = 1,2,3,4, 



(1.13) 



for ui > M2 > U3 > W4. In Figure 1, we plot the self-similar solution of the Whitham 
equations (1.7) with g ~ I for the NLS equation, and the corresponding periodic 
oscillatory solution (1.8) for the initial data (1.11) with a = 4, 6 = and c = —1. The 
oscillations describe a dispersive shock of the NLS equation under a small dispersion. 
Note here that the oscillations have a uniform structure, which is due to an almost 
linear profile of the Whitham solution U3 . This will be seen to be in sharp contrast 
to the case of the mKdV equation, which we will discuss later (cf. Figure 1.2). 



4 
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Fig. 1.1. Self-Similar solution of the NLS-Whitham equation (1.7) of g = 1 and the corre- 
sponding oscillatory solution (1.8) of the NLS equation with e = 0.15. The dark line in the middle 
of the oscillations is the weak limit p{x,t) given by (1.10). The initial data is given by (1.11) with 
a = 4, b = and c = — 1. 



The figures in this paper all have the same form: On the left hand side is a plot 
of the solution of the Whitham equations as a function of the self-similarity variable 
x/t, which is exact, other than a numerical method used to implement the inverse 
function theorem. On the right hand side is the oscillatory solution given by (1.8) 
(respectively (1.21) for mKdV) a± t ~\, while the dark plot is the weak limit (1.10) 
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of the oscillatory solution at t = 1, both plots on the right are also exact. In the 

first two figures wc dcniark the region where the Whithani equations with g = I 
govern the solution, place a dashed-dotted line where the behavior of the oscillatory 
solution changes, and label the four functions ui > U2 > U3 > U4. The demarcation 
and labeling are similar in the other figures and we will leave them off for brevity. 
Although we do not include any numerical simulations, we would like to mention that 
E. Overman showed us his numerical simulation of the NLS and mKdV equations 
which captures the features of the oscillatory solutions plotted in this paper. 

The defocusing NLS equation is just the first member of the defocusing NLS 
hierarchy; the second is the (defocusing) complex modified KdV (mKdV) equation 



dt 2 



2^ 

dx 



= 



We again use V'(a:, t; e) = A{x, t; e) exp (^^/^ ^l^M) j g^j^^^ notation p{x, t; e) = A'^{x, t; e), 
v{x, t; e) = dS{x, t; e)/dx to obtain the conservation form of the mKdV equation 



4 dx^ 

six,t,e) ^ notation p{x, t; e) 



(1.14) 



dp 9/89 3 . 
di + d-x (4^ + -4^' 

d. , d 

dt^'^'^ + dx 



dx 



-p v+-pV^ 



,3/4 

d_ 
dx 



dx"^' 



.1/4 



d'^ 3 ■ 



(L15) 



where 



R = 



3v f dp 



2p \dx 



dv dp 
dx dx 



d^p 

dx/ 



The mass density p and momentum density pv for the mKdV equation also have 
weak limits as e ^ [6]. As in the NLS case, the weak limits satisfy 



d_ d_ 

di^^ di 



{pv) 



d_ 
dx 



3 2 
p V - 



3 3 
-pv 



(1.16) 



, 



dV''"' ' dx\2'' " ' A' 
until the solution of (1.16) forms a shock. One can rewrite equations (1.16) as 



d_ (a 
dt \I3 



3 fba^ 



2a(3 + 0^ \ 



, 



(1.17) 



where the Riemann invariants a and (3 are again given by formula (1.5). 

Let us again consider the simplest case a{x, 0) = a, where a is constant, to see 
how the solution of system (1.17) develops a shock. In this case, the system reduces 
to a single equation 



9/3 



d(3 



dt ' 8(«^ + 2a/3 + 5/3^)^=0. 

As in the NLS case, we consider the initial function given by /3(.x, 0) = for a; < 
and /3(x, 0) = c for a; > 0. We recall that, in the NLS case, the zero phase solution of 
(1.4) develops a shock if and only if 6 > c. However, the solution in the mKdV case 
develops a shock for 6 > c if and only if a + 56 > 0. In addition, if 6 < c, the solution 
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in the mKdV case develop a shock if and only if a + 56 < 0. These differences between 
the niKdV and NLS cases are due to the weak hypcrbolicity of the system (1.17) (note 
that for the eigenspeed A = + 2a(3 + 5/3^) for /?, we have d\/d(} = |(a + 5/3) 
which can change sign). As will be shown below, this leads to an additional structure 
in the dispersive shock for the niKdV case. 

As in the case of the NLS equation, immediately after the shock formation in the 
solution of (1.16), the weak limits are described by the mKdV-Whitham equations 

^+M9,i(Mi,...,W2g+2)^ =0 , i = l,2,3...,25 + 2 , (1.18) 

where fig j's can also be expressed in terms of complete hyperelliptic integrals of genus 
9 [6]. 

In this paper, we study the solution of the Whitham equations (1.18) with g =1 
when the initial mass density p(x, 0) and momentum density p{x, 0)f (x, 0) are step 
functions. In view of (1.5), this amounts to requiring a and /? of system (1.17) to 
have step-like initial data. We are interested in the following two cases: 

(i) a{x, 0) is a constant and 

{h X < 
^>0' a>b ,a> c ,by^ c , (1.19) 

(ii) (3{x, 0) is a constant and 

a(x,0) = |^' !^ ^ Q 7 /3(x, 0) = a, b> a ,c> a ,b^ c . (1.20) 

In the case of the NLS equation, the genuine nonlinearity of the single phase 
Whitham equations (see (1.13)) warrants that the solution is found by the implicit 
function theorem. However the mKdV-Whitham equations (1.18), in general, are not 
genuinely nonlinear, that is, a property like (1.13) is not available (see Lemma 2.1 
below). Our construction of solutions of the Whitham equation (1.18) with g = 1 
makes use of the non-strict hypcrbolicity of the equations. For the NLS case, it 
has been known in [6, 7] that the Whitham equations (1.7) with <? = 1 are strictly 
hyperbolic, that is, 

for ui > 1*2 > M3 > M4. For the mKdV-Whitham equations (1.18) with g = 1, the 
eigenspeeds ^i,i(wi, U2, W3, W4) niay coalesce in the region Mi > U2 > W3 > W4 (we will 
discuss the details in Section 2). 

Let us now describe one of our main results (see Theorem 3.1) for the single phase 
mKdV- Whitham equations with step- like initial function (1.19) for a = 4, b = 1 and 
c = —1. In this case, the space time is divided into four regions (see Figure 1.2) 
instead of three in the case of the NLS equation (cf. Figure 1) 

(1) I < ci , (2) ci < I < C2 , (3) C2 < I < cs , (4) I > C3 , 

where Ci, C2 and C3 are some constants. In the first and fourth regions, the solution 
of the 2x2 system (1.17) governs the evolution: 

(1) for x/t < ci. 



a{x,t) = 4, P{x,t) = 1, 
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(4) for x/t > C2, 

a.{x,t)=A, I3{x,t) = -1. 

The Whitham solution of the 4x4 system (1.18) with g = 1 lives in the second and 
third regions; 

(2) for ci < x/t < C2, 

X X 
Ui{x,t)=A, U2{x,t) = l, - = /Xi,3(4, 1,M3,M4), - = /Zi,4(4, 1, M3, W4) , 

(3) for C2 < x/t < C3, 

X 

Ui{x,t)=A, U2{x,t) = l, - = /ii,3(4,l,U3,-l) , U4{x,t) = -1. 

Note that, in the second region, we have 

/il,3(4, 1, U3, U4) = /(Xl,4(4, 1, ^3,^4) 

on a curve in the region — 1 < M4 < U3 < 4. This implies the non-strict hyperbolicity 
of the mKdV- Whitham equations (1.18) for 5' = 1. 

It is again possible to use the method of [1, 16] to show that the solution of the 
mKdV equation (1.14) can be approximately described, in the single phase regime, 
by the periodic solution of the mKdV when e is small. The periodic solution has the 
same form as (1.8) of the NLS, i.e., 

p{x,t;e) = P3 + {p2 - Pa) sia^WPi - P3 t; e), s) . (1.21) 

However, 6{x,t;e) is now given by ^ = (a; — V2t)/e with the velocity V2 (see e.g. [9]) 

where cti = X]j=i '^j ^^^^ "'2 = Z^i<j UiUj are the elementary symmetric functions of 
degree one and two, respectively. The functions pi , p2 and p3 are also given by formula 
(1.9). If Ml, U2, W3 and U4 are constants, formula (1.21) gives the periodic solution 
of the mKdV equation. To describe the solution p{x,t;e) of the mKdV equation 
(1.14), the quantities Ui, U2, W3 and U4 must satisfy the single phase mKdV- Whitham 
equations (1.18) for g = 1. The weak limit of p{x,t; e) of the mKdV equation is also 
given by formula (1.10). 

In Figure 1.2, we plot the self-similar solution of the Whitham equations (1.18) for 
g = 1 and the corresponding periodic oscillatory solution (1.21). We note here that 
the pattern of the oscillation in this case has two distinct structures: one corresponds 
to the region (2), ci < x/t < C2, and the other corresponds to the region (3), C2 < 
x/t < C3, We also note that the weak limit p{x,t) is not smooth at the boundary 
point x/t = C2 ~ 4.63. 

As we will show below, for other values of a, b and c, the solutions of (1.17) and 
(1.18) with g = 1 will be seen to be quite different from the above. 

The Whitham equations (1.18) for the mKdV equation are analogous to the 
Whitham equations for the higher members of the KdV hierarchy [11, 12]. For step 
like initial data, the single phase Whitham solutions for the higher order KdV are also 
constructed using the non-strict hyperbolicity of the equations. In the case of strictly 
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Fig. 1.2. Self-Similar solution of the mKdV-Whitham equation (1.18) with g = 1 and the 
corresponding oscillatory solution (1.21) of the mKdV equation with e = 0.012. The initial data is 
given by (1.19) with a = 4, 6 = and c = —1. There are two distinct structures in the oscillations 
and they are separated by x/t ^ 4.63 (cf. Figure 1). 



hyperbolic Whitham equations for the KdV, the oscillations (dispersive shock) have 
uniform structure. However, in the case of non-strict hyperbolic Whitham equations 
for the higher order KdV, an additional structure has been found in the dispersive 
shocks. This new structure is similar to the one found here in the dispersive shocks 
of the mKdV equation. 

The organization of the paper is as follows. In Section 2, we will study the 
eigenspeeds, fig^i, ftg,2, ^J■g,3 and ftg^4 of the Whitham equations (1.18) for g — 1. 
In Section 3, we will construct the self-similar solutions of the single phase Whitham 
equations for the initial function (1.19) with a > b > c. In Section 4, we will construct 
the self-similar solution of the Whitham equations for the initial function (1.19) with 
a > c > b. In Section 5, we will briefly discuss how to handle the other step- like initial 
data (1.20). 

2. The Whitham Equations. In this section we define the eigenspeeds Xg/s 
and /ig,i's of the Whitham equations (1.7) and (1.18) with g = I for the NLS and 
the mKdV equations. For simplicity, we suppress the subscript g = 1 in the notation 
Xg/s and ^J,g/s in the rest of the paper. 

We first introduce the polynomials of ^ for n = 0, 1, 2, . . . [2, 4, 10]: 

Pn{^,Ui,U2,U3,U4) = + an,lC + ■ ■ ■ + an.n+2 , (2.1) 

where the coefficients, a„^i, a„_2, • ■ ■ , (in,n+2 are uniquely determined by the two con- 
ditions 

P„{£.,Ui,U2,U3,Ui) /o^C-2n f 1 Id 

^ +0{S, ) for large \i\ 

V (C - ■"l)(C - U2){i - U3)(C - "4) 



and 



P„(^,Wl,U2,U3,U4) 



«2 V ("1 - - ■"2)(C - "3)(^ - Ui) 



The coefficients of P„ can be expressed in terms of complete elliptic integrals. 

The eigenspeeds of the Whitham equations (1.7) with g = 1 for the NLS equation 
are defined in terms of Pq and Pi of (2.1) [4, 6, 10], 

w N s Pl{Ui,Ul,U2,U3,U4) . -, „ „ . 

Ai(Ul,U2,M3,W4) = 8 — r, 1=1,2,3,4, 

Po{Ui,Ui,U2,U3,Ui) 
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which give 

Xi{Ui,U2,U3,U4) = 2 ai{Ui,U2,U3,U4) - — — ^ . (2.2) 

Here ai := ^^-lUj, and M2, U3, U4) is given by a complete elliptic integral [13] 

/""' dn 
I{uuU2,U3,U4)= . (2.3) 

Ju2 \J {ui - ri)[ri - U2 ){r] - U3)(r] - U4) 

The function / can be rewritten as a contour integral. Hence, 

dl dl 

2("^-"j)7r^= ' «,i = 1,2,3,4 (2.4) 

OUiOUj OUi OUj 

since the integrand satisfies the same equations for each -q ^ Ui. This contour integral 
connection also allows us to give another formulation of / 

I{uuU2,U3,U4) = . (2.5) 

Jui \/{ui-ri) (U2 - V) [Us -ri){r]-U4) 

It follows from (2.2), (2.3) and (2.5) that 

A4 - 2c7i < A3 - 2c7i < < A2 - 2cti < Ai - 2ai (2.6) 

for 7i4 < M3 < 1*2 < Ui. This implies the strict hyperbolicity of the NLS-Whitham 
equation (1.7) for = 1. 

The eigenspeeds Aj's have the following values [13]: At U3 = U4, we have 

Ai = 6ui + 2u2, 

A2 = 2ui+6u2, (2.7) 

A3 = A4 = 2(ni +U2+ 2u4) - '^x:tx:'\ 

and at W2 = U3, 

Ai = 6ui + 2u4, 

A2 = A3 = 2ui + 4u3 + 2u4, (2.8) 
A4 = 2ui + 6U4. 

Notice that the cigcnspccd A2 = A3 at U2 = U3 is the same as the velocity of the 
periodic solution (1.8), i.e. Vx = 2u\ = 2{u\ + 2u3 + U4). 

The eigenspeeds of the mKdV-Whitham equations (1.18) with g = \ are [6] 

, ^ o -f2(Wi,Ml,«2,W3,U4) ■ i o Q /I n^ 

/Xi(ui,U2,W3,W4) = 3 ^-^ r- , 1=1,2,3,4. (2.9) 

Po(Mi, Wl,'"2,W3,U4) 

They can be expressed in terms of Ai, A2, A3 and A4 of the NLS-Whitham equations 
(1.7) with ,g = 1. 

Lemma 2.1. The eigenspeeds iJi{ui,U2,U3,U4) 's of (2.9) can he expressed in the 
form 

= i(Ai-2ai)^+g , z= 1,2,3,4, (2.10) 
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where a\ = % ^"^^ 9 — <l{ui,U2,U3,U4) is the solution of the boundary value 

problem of the Euler-Poisson-Darboux equations 



d^q dq dq 
duiduj dui duj ' 
q{u, u, u) = 3u^. 



2K-%)l.r^ = ^-l^ ' i,j = 1,2,3,4, (2.11) 



Also the /ij 's satisfy the over-determined systems 

1 dfXi 1 dXi 

jjLi — jjLj duj Xi — Xj duj 



i^j . (2.12) 



We omit the proof since it is very similar to the proof of an analogous result for 
the KdV hierarchy [14]. 

The boundary value problem (2.11) has a unique solution. The solution is a 
symmetric quadratic function of Mi, U2, Us and U4 

^^^2 , (2.13) 

where a2 = J2i>j ^i"^] is the elementary symmetric polynomial of degree two. Notice 
that q gives the velocity of the periodic solution (1.21) for the mKdV equation, i.e. 
V2 = q. 

For NLS, A,'s satisfy [13] 

5A4 ^ 3A3-A4 ^5A3 
dui 2us — U4, dus 

for U4 < U3 < U2 < Ui. Similar results also hold for the mKdV-Whitham equations 
(1.18) with g = l. 



Lemma 2.2. 



p>l^^ ^f^>0, (2.15) 

dn4 ^3 fi4 .f dq 
au4 2 M3 — M4 0U4 



for U4 < U3 < U2 < ui. 

Proof We use (2.10) and (2.14) to obtain 

a/X3 ^ 1 aAa _^ 1 _2^)^ 
dus 2du3du3 2 ^ ^ 9u'^ 



and 



^ 3 X3-X4 dq , 1^^ o^^^^9 ro^7^ 
>7 1- -(A4 - 2ai)— ^ , (2.17) 



1 / X \ \ dq I ^ \( dq dq 

M3-M4=2(A3-A4)^ + 2(A4-2a0(^^-^ 

1 . ^ dq „ X/ N d'^q 

= o ('^^S - A4) ^ h (A4 - 2(7i)(U3 - U4)- — 

2 du3 OU3OU4 

2, , /3 A3- A4 dq 3,, „ ^ oi2g , 
= -{U3 -U4)[- + - A4 - 2ai — 2- , 2.18 

3 V4 W2 — U3 0U3 2 OU3OU4 ' 
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where we have used equation (2.11) 



dq 



— = 2(U3 - ^4) ^ ^ ■ 



It follows from formula (2.13) for q that 



which, along with with (2.17) and (2.18), proves (2.15). Inequahty (2.16) is proved in 
the same way. □ 

The following calculations are useful in the subsequent sections. Using formula 
(2.10) for 113 and /Z4 and formulae (2.2) for A3 and A4, we obtain 



Ms - M4 = 



I 


dq 


dl 


dq 


dl 




du4 






du4 



{duj){duj) 

2J(u3 — U4) 



dq dl 
dui du3 



M 



dl . _ , dq 
dui du3 



dq dl 
dui dui 



(2.19) 



where 



M = 



dq d^I 



d^q dl 



dui du^dui duT,du4 du^ 

Here we have used equations (2.4) for / and equations (2.11) for q in equality (2.19). 
Since q of (2.13) is quadratic, we obtain 



dM _ dq d^I 
dus dui du^dui 



(2.20) 



We note that another expression for M is 

dq d^I 



M = 



d^q dl 



Hence, we get 



du3 du^dui dusdui dus 



dM _ dq d^I 
du4 dus dusdul 



(2.21) 



We next evaluate M(ui,it2,«3,«4) when U3 = U4. Using the integral formula 
(2.3) for the function I and applying the change of variable r] = (ui — U2)i' + U2, we 
obtain 



M 



du4 



1{u2-U4Yk {l+^vf^^Al^) 
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The two integrals can be evaluated exactly as 

dp tt{2 + 7) 



dp _ 7r(8 + 87 + 87^) 

Jo il + -fp)^y/p{l-p) ~ 8(1 + 7)5 ' Jo il + -/iy)^^ypil-p) 2(1 + 7)1 
for 7 > —1. We finally get 

TrU{ui,U2,U4) 



M 
where 



128[(m2 — U4){ui — W4)]- 



(2.22) 



U{ui,U2,0 = [8(w2 - + 8(u2 - e)(Mi - U2) + 3(?ii - U2f]{ui +U2+ 4^ 
- 8{U2 - 0[2(W2 - 0^ + 3{U2 - 0(wi - U2) 

+{ui-U2f]. (2.23) 
Similar to (2.19) for jj^ and /Z4, we have 

where 

dq d'^I d'^q dl 

N = — . 

du2 du2du3 du2du3 du2 

Since q of (2.13) is quadratic, wc obtain 

dN dq d^I 



dus du2 du2du\ 
Finally, we use (2.7) and (2.10) to calculate 



(2.25) 



1 1 

(M2 - Ms) = 2(U1 + U2 + 2U4)]^ oI'^S - 2(ui + U2 + '^Ui)\- 

U3=M4 ^ OU2 ^ 

= or ^ V{u„U2,u,) , (2.26) 
2{Ui + U2 — 2U4) 

where 

F(ui, ^2, ^4) = 3uf + 3u| — 12^4 + 6uiU2 + 6U1U4 — 6U2W4 . (2.27) 

3. Self-Similar Solutions. In this section, we construct self-similar solutions of 

the Whitham equations (1-18) with g = 1 for the initial function (1.19) with a > b > c. 
The case with a > c > b will be studied in next section. The solution of the zero 
phase Whitham equations (1.17) does not develop a shock when a + 56 < 0. We are 
therefore only interested in the c;asc a + 5b > 0. 

We first study the ^-zero of the cubic polynomial equation 

[/(a,6,0=0, (3.1) 



where U is given by (2.23). It is easy to prove that for each pair of a and b satisfying 
a> b and o + 56 > 0, U{a, 6, ^) =0 has only one simple real root. Denoting this zero 
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by ^(a, 6), we then deduce that U{a,b,^) is positive for ^ > ^{a,b) and negative for 
^ < ^(a, 6). Since U{a, 6, — (a + 6)/4) < in view of (2.23), we must have 

^(a,b)>-^. (3.2) 

For initial function (1.19) with a> b> c and a + 56 > 0, we now classify the resulting 
Whitham solutions into four types: 

(I) ^{a, b) < c with any a> b > c 
(II) lla, b)> c with a + 56 > 3(6 - c) > 

(III) l{a, b)> c with a + 56 = 3(6 - c) > 

(IV) l{a, b)>c with < a + 56 < 3(6 - c) 
We will study the second type (II) first. 

3.1. Type II. Here we consider the step initial function (1.19) satisfying ^(o, 6) > 
c and a + 56 > 3(6 - c) > 0. 

Theorem 3.1. (see Figure 1.2.) For the step-like initial data (1.19) with a > 

b > c, a + 56 > 3(6 — c) and £,{a,b) > c, the solution (a,/?) of the zero phase Whitham 
equations (1.17) and the solution (wi, U2, M3, M4) of the single phase Whitham equations 
(1.18) with g = 1 are given as follows: 

(1) For x/t < M3(a, 6, ^{a, b),^{a, b)), 

a = a, l3 = b. (3.3) 

(2) For /L<3(a,6,C(a,6),C(a,6)) < x/t < fi3{a,b,u** ,c), 

ui=a, U2 = b, - = i^s{a,b,U3,U4) , - = /i4(a, 6, us, U4) , (3.4) 

where u** is the unique solution U3 of H3{a,b,U3,c) = /i4(a, 6, U3, c) in the 

interval c < < b. 

(3) For iJ,3{a, 6, u**,c) < x/t < ^^{a, 6, 6, c), 

X 

ui = a , U2 = b , - = /X3(a, 6, Us, c) , U4 = c. (3.5) 

(4) For x/t > fi3{a, 6, 6, c), 

a = a , P = c. (3-6) 

The boundaries x/t = /X3(a, 6, ^(o, 6), ^(a, 6)) and x/t = /X3(o, 6, 6, c) are called 
the trailing and leading edges, respectively. They separate the solutions of the single 
phase Whitham equations (1.18) with (/ = 1 and the zero phase Whitham equations 
(1.17). The single phase Whitham solution matches the zero phase Whitham solution 



in the following fashion (see Figure 1.): 

(ui,M2) = the solution (a,/3) of (1.17) defined outside the region , (3.7) 

U3 = Ui , (3.8) 

at the trailing edge; 

{ui,Ui) = the solution (a, /3) of (1.17) defined outside the region , (3.9) 

U2 = U3 , (3.10) 



14 



Y. KODAMA, V.U. PIERCE AND F.-R. TIAN 



at the leading edge. 

The proof of Theorem 3.1 is based on a series of lemmas: We first show that the 
solutions defined by formulae (3.4) and (3.5) indeed satisfy the Whitham equations 

(1.18) for g = 1 [2, 15]. 
Lemma 3.2. 

(1) The functions Ui, U2, W3 and U4 determined by equations (3.4) give a solution 
of the Whitham equations (1-18) with g = 1 as long as U3 and u^ can he solved 
from, (3.4) as functions of x andt. 

(2) The functions ui, U2, M3 and U4 determined by equations (3.5) give a solution 
of the Whitham equations (1-18) with g =1 as long as U3 can be solved from 
(3.5) as a function of x andt. 

Proof. (1) Ui and U2 obviously satisfy the first two equations of (1.18) for g = I. 
To verify the third and fourth equations, we observe that 

p = p=0 (3.11) 

OU4 OUs 

on the solution of (3.4). To see this, we use (2.12) to calculate 

« 9 A3 

The second part of (3.11) can be shown in the same way. We then calculate the partial 
derivatives of the third equation of (3.4) with respect to x and t 

1 = ^ tUSx , = tUst + M3 , 

which give the third equation of (1.18) with g = 1. The fourth equation of (1.18) with 
g = 1 can be verified in the same way. 

(2) The second part of Lemma 3.2 can easily be proved. □ 

We now determine the trailing edge. Eliminating x and t from the last two 
equations of (3.4) yields 

/i3(o, b, U3, U4) - 114(0, b, U3, U4) = . (3-12) 

Since it degenerates at U3 = U4, we replace (3.12) by 

rp/^h. . ^ M3(a, b, Us, U4) - M4(a, b, U3, U4) 
[a, 0, U3,U4) := = . (3.13) 

U3 - U4 

Therefore, at the trailing edge where U3 = U4, equation (3.13), in view of formulae 

(2.19) and (2.22), reduces to 

C/(a,6,U4) = . (3.14) 

Noting that ^(a, h) is the unique solution of (3.1), we then deduce that U4 = ^(o, b). 

Lemma 3.3. Equation (3.13) has a unique solution satisfying U3 = U4. The 
solution is U3 = U4 = ^{a,b). The rest of equations (3.4) at the trailing edge are 
u\ = a, U2 = b and x/t = H3{a, b,^{a, b),^{a, b)). 

Having located the trailing edge, we now solve equations (3.4) in the neighborhood 
of the trailing edge. We first consider equation (3.13). We use (2.19) to write F of 
(3.13) as 

27 

F{a,b,U3,U4) = M{a,b,U3,U4) ■ 
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We note that at the traUing edge U3 = U4 = ^{a, b), we have M{a, b, ^{a, b),^{a, b)) = 
because of (2.22) and (3.14). We then use (2.20) and (2.21) to differentiate F at the 
traihng edge 

dF{a, b, e(a, &), e(a, b)) dF{a, b, ^{a, b),^{a, b)) 



dus du. 



■4 



[a + 6 + 4^(a,6)] ^^-2^>0, 



where we have used the expression (2.13) for q in the last equation and (3.2) in the 
inequality. These show that equation (3.13) or equivalently (3.12) can be inverted to 
give W4 as a decreasing function of U3 

Ui = A{us) (3.15) 

in a neighborhood of 1/3 = 1*4 = ^(a, b). 

We now extend the solution ^(^3) of equation (3.12) in the region c < W4 < 
^(a, 6) < Us < 6 as far as possible. We first claim that 

dq{a,b,U3,U4) ^ ^ dq{a,b,U3,U4) ^ „ . . 

on the extension. To see this, we first observe that inequalities (3.16) are true at 
the trailing edge tt3 = ^4 = ^(a, 6). This follows from (2.13) and (3.2). Therefore, 
inequalities (3.16) hold in a neighborhood of the traihng edge. To prove that (3.16) 
remains true on the extension, we use formula (2.10) for 113 and ^4 to rewrite equation 
(3.12) as 

-[A3 -2(a + 6 + W3 + M4)]7r^ = -[X4-2{a + b + U3 + U4)]^ . 
z 0U3 z 0U4 

Since the two terms in the two parentheses are both negative in view of (2.6) and 
since ^ — ^ = (^3 — U4)/4 > in view of (2.13), neither ^ nor can vanish 
on the extension. This proves inequalities (3.16). 
We deduce from Lemma 2.2 that 

P>0, p<0 (3.17) 
0U3 0U4 

on the solution of (3.12). Because of (3.11) and (3.17), solution (3.15) of equation 
(3.12) can be extended as long as c < U4 < ^(a, b) < U3 < b. 

There are two possibilities: (1) U3 touches b before or simultaneously as U4 reaches 
c and (2) U4 touches c before U3 reaches b. It follows from (2.8), (2.10) and (2.13) 
that 

IJ.3{a,b,b,U4) — iJ,4{a,b,b,U4) = ^ {b — U4){a + 2b + 3u4) > for c < U4 < 6 , (3.18) 

where we have used a + 26 + 3c > in the inequality. This shows that (1) is unattain- 
able. Hence, U4 will touch c before U3 reaches b. When this happens, equation (3.12) 

becomes 



H3{a,b,U3,c) = iJ,4{a,b,U3,c) 



(3.19) 
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Lemma 3.4. Equation (3.19) has a simple zero in the interval c < uz < b, 
counting multiplicities. Denoting the zero by u**, then /U3(a, 6, M3, c) — jj,4,{a,b,U3,c) 
is positive for U3 > u** and negative for U3 < u** . 

Proof. We use (2.19) and (2.20) to prove the lemma. In both formulae, du^I, 
du^I and d"^^.^^! are all positive functions. By (2.20), 

dM{a,b,ua,c) (a + 6 + U3 + 3c) 

hi = A a 2a for c < W3 < 6 . (3.20) 

We claim that 

M{a, b, M3, c) < when U3 = c and M{a, b,us,c) > for us near b . 

The second inequality follows from (2.19) and (3.18). The first inequality can be 
deduced from formula (2.22) 

^' ^) = ^^^^' ^'"^ < for c < b) . 
128[(o — c)(a — c)J 2 

Therefore, Af (a, 6, 713, c) has a zero in the interval c < U3 < 6. The uniqueness of the 
zero follows from (3.20) in that M{a, b, U3, c) increases or changes from decreasing to 
increasing as uz increases. This zero is exactly u** and the rest of the theorem can 
be proved easily. □ 

Having solved equation (3.12) for U4 as a decreasing function of U3 for c < U4 < 
^(0,6) < U3 < 6, we turn to equations (3.4). Because of (3.11) and (3.17), the third 
equation of (3.4) gives M3 as an increasing function of x/t, for /i3(a, b, ^(a, 6), ^(a, b)) < 
x/t < ijL3{a,b,u** ,c). Consequently, U4 is a decreasing function of x/t in the same 
interval. 

Lemma 3.5. The last two equations of (3.4) can be inverted to give and U4 
as increasing and decreasing functions, respectively, of the self- similarity variable x/t 
in the interval iJ,s{a,b,^{a,b),^{a,b)) < x/t < iJLz{a,b,u** ,c), where u** is given in 
Lemma 3.4. 

We now turn to equations (3.5). We want to solve the third equation when 
x/t > /i3 (a, &,«**, c) or equivalently when M3 > u**. According to Lemma 3.4, 

113(0, b, U3, c) — /i4(a, 6, U3, c) > for u** < U3 < b. In view of (3.16), du^q{a, 6, U3, c) = 
(a + 6+3u3 + c)/4 is positive at U3 = u** and hence, it remains positive for U3 > u**. 
By (2.15), we have 

dn3{a,b,U3,c) ^ Q 

dU3 

Hence, the third equation of (3.5) can be solved for U3 as an increasing function of 
x/t as long as M** < M3 < b. When U3 reaches 6, we have x/t = /U3(a, 6, 6, c). We have 
therefore proved the following result. 

Lemma 3.6. The third equation of (3.5) can be inverted to give U3 as an increasing 
function of x/t in the interval 113(0, b, u** , c) < x/t < 113(0, b,b,c). 

We are ready to conclude the proof of Theorem 3.1. The solutions (3.3) and (3.6) 
are obvious. According to Lemma 3.5, the last two equations of (3.4) determine U3 
and U4 as functions of x/t in the region 113(0, b, ^(a, b), ^(a, b)) < x/t < /i3(a, b, u** , c). 
By the first part of Lemma 3.2, the resulting ui, U2, U3 and U4 satisfy the Whitham 
equations (1.18) with g — I. Furthermore, the boundary conditions (3.7) and (3.8) 
are satisfied at the trailing edge x = ii3(a, b, ^(a, b),^(a, b)). 
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Similarly, by Lemma 3.6, the third equation of (3.5) determines as a function of 
x/t in the region /i3(a, b, m**, c) < x/t < ^i^i^a, b, b, c). It then follows from the second 
part of Lemma 3.2 that ui, U2, U3 and 7/4 of (3.5) satisfy the Whitham equations 
(1.18) for g = I- They also satisfy the boundary conditions (3.9) and (3.10) at the 
leading edge x/t = ^3{a, b, b, c). We have therefore completed the proof of Theorem 
3.1. 

A graph of the Whitham solution (ui, U2, U3, U4) is given in Figure 1.2. It is 
obtained by plotting the exact solutions of (3.4) and (3.5). 

3.2. Type I. Here we consider the initial function (1.19) satisfying £,{a,b) < c 
with b > c and a + 56 > 0. 




Fig. 3.1. Self-Similar solution of the Whitham equations (1.18) with g = 1 and the correspond- 
ing oscillatory solution (1.21) of the mKdV equation with e = 0.07. The initial data is given by 
(1.19) with a = 7 , b = and c = —1 of type I. 



We will only present our proofs briefly, since they are, more or less, similar to 
those in Section 3.1. The main feature of this case is that the ^-zero point does not 
appear in the solution u^, and the Whitham equations (1.18) with g = I are strictly 
hyperbolic on the solution. 

Theorem 3.7. (see Figure 3.1.) For the step-like initial data (1.19) with a > 
b > c, a + 56 > and ^(a, 6) < c, the solution of the Whitham equations (1.18) with 
g — 1 is given by 

X 

Ml = a , U2 = 6 , - = ^3(0, 6,M3,C) , -U4 = c 

for fi3{a,b,c,c) < x/t < ii3(a,b,b,c). Outside this interval, the solution of (1.17) is 
given by 

X 

a = a, j3 = b for - < H3{a,b,c,c) , 

and 

X 

a = a , P = c for — > ^3{a,b,b,c) . 



Proof. It suffices to show that fj,3(a,b,U3, c) is an increasing function of U3 for 
c < U3 < 6. Substituting (2.13) for q into (2.20) yields 

dM(a,b,U3,c) K , d^I I, , d"^! 



18 



Y. KODAMA, V.U. PIERCE AND F.-R. TIAN 



for c < U3 < 6, where we have used ^(a, 6) < c in the first inequality and (3.2) in 
the second one. We now use formula (2.22) to calculate the value of M(a, b, U3, c) at 

U3 = c 

TrU{a, b, c) 



M{a,b,c, c) = 



> for C(a, b) <c 



128[(6-c)(a-c)]5 

because U{a,b,£,) > for ^ > S,{a,b). Therefore, M{a,b,U3,c) > for c < 1*3 < fe. It 
then follows from (2.19) that /Lt3(a, b, U3, c) — fj,4{a, b, U3, c) > 0. Since J^(a, b, U3, c) = 
(o+6+3m3 + c)/4 > (a + 6+3^(a, 5))/4 > because of (3.2), we conclude from Lemma 
2.2 that 

d//3(a, &,M3,c) 



du-. 



> 



for c < M3 < 6. □ 

3.3. Type III. Here we consider the step initial function (1.19) satisfying ^(a, b) > 
c with a + 5& = 3(& - c) > 0. 



3 
2 



-1 



2.4 



P(x,1) 



2.6 

x/t 



2.8 




Fig. 3.2. Self-Similar solution of the Whitham equations (1.18) with g = 1 and the correspond- 
ing oscillatory solution (1.21) of the mKdV equation with e = 0.007. The initial data is given by 
(1.19) with a = 3, b = 0, c = -1 0/ type III. 

Theorem 3.8. (see Figure 3.2.) For the step-like initial data (1.19) with a > 
h > c, ^(a, b) > c and a + 56 = 3(6 — c), the solution of the g = I Whitham equations 
(1.18) with g ~ 1 is given by 

X X 

ui = a , W2 = 6 , y = /i3(a, 6, U3,U4) , - = ^4(0, 6, U3, M4) , 

for /X3(a, 6, ^(a, 6), ^(a, 6)) < x/t < fi3(a,b,b,c). Outside the region, the solution of 
equations (1.17) is given by 



and 



/3 = 6, for J < fi3{a,b,^{a,b),^{a,b)) . 



a = a, (3 = c, for - > ^i3{a,b,b,c) . 



Proof. It suffices to show that M3 and U4 of /X2(a, 6, 1*3, U4) — /i3(a, 6, U3, U4) = 
reaches 6 and c, respectively, simultaneously. To see this, we deduce from the 
calculation (3.18) that 

//3(a, 6, 6, U4) — ^i{a, 6, 6, M4) = -(6 — M4)(a + 26 + 3^4) , (3.21) 



vanishes at W4 = (—a — 26)/3 = c. □ 
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3.4. Type IV. Here we consider the step initial function (1.19) satisfying ^(a, 6) > 
c with < a + 56 < 3(6 - c) . 




p(x,1) 




Fig. 3.3. Self-Similar solution of the Whitham equations (1.18) with g = 1 and the correspond- 
ing oscillatory solution (1.21) of the mKdV equation with e = 0.006. The initial data is given by 
(1.19) with a = 2.75, b = and c = —1 of type IV. The solution in the region 2.52 < x/t < 2.65 
represents a rarefaction wave. 

Theorem 3.9. (see Figure 3.3.) For the step-like initial data (1.19) with a > 
b > c. ^(a, b) > c and < a + 56 < 3(6 — c), the solution of the Whitham equations 
(1.14) is given by 



ui = a , U2 



- = ^J■2{a, 6, U3, Ui) , x = ^i3{a, 6, M3, U4) t 



for fj,3{a,b,^{a,b),^{a,b)) < x/t < /X3(a, 6, 6, — (a + 26)/3). Outside the region, the 
solution of equation (1.17) is divided into the three regions: 

(1) For x/t < ti:i{a,b,i{a,b),i{a,b)), 

a = a , (3 ~ b . 

(2) For fi3{a,b,b,-{a + 2b)/3) < f < | {a^ + 2ac + 5c^) , 

1 



a = a, l3 
(3) For x/t> ^ (a'^ + 2ac + 5c^), 



8 X 
15 t 



[3 = c. 
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Proof. By the calculation (3.21), when U3 of ii3(a,b,U3,Ui) — /i4(a, 6, U3, U4) = 
touches 6, the corresponding U4 reaches — (a + 26)/3, which is above c. Hence, 
equations 

X X 

- = fL3{a,b,U3,U4) , - = /i4(a, 6, ■ti3,"4) 

can be inverted to give U3 and U4 as functions oix/t in the region /X3(a, 6, ^(a, 6), ^(a, 6)) < 
x/t < ^^3(0, 6, 6, — (a + 26)/3). In the region (2), the equations (1.17) has a rarefaction 
wave solution. □ 

4. More Self-Similar Solutions. In this section, we construct self-similar so- 
lutions of the g = I Whitham equations (1.18) for the initial function (1.19) with 
a > c > b. The solution of equations (1.17) does not develop a shock for a -I- 56 > 0. 
We are therefore only interested in the case a -I- 56 < 0. We classify the resulting 
Whitham solution into four types: 
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(V) a + 56 < -4(c- &) < 
(VI) > a + 56 > -4(c - b) with V{a, c, 6) < 
(VII) > a + 56 > -4(c - b) with V{a, c, 6) = 
(VIII) > a + 56 > -4(c - 6) with Via, c, 6) > 
where y is a quadratic polynomial given by (2.27). 

4.1. Type V. Here we consider the step initial function (1.19) satisfying a+56 < 
-4(c- 6) < 0. 



0.5 


-0.5 
-1 




1 2 3 1 2 3 

x/t X 



Fig. 4.1. Self-Similar solution of the Whitham equations (1,18) with 3 = 1 and the correspond- 
ing oscillatory solution (1.21) of the mKdV equation with e = 0.018. The initial data is given by 
(1.19) with a = 1/2, b = -I and c = of type V. 



Theorem 4.1. (see Figure ^.l.) 
c > b, a + 56 < — 4(c — 6), the solution 
given by 

ui = a , U2 — c , 

for fi-i{a,c,c,b) < x/t < ^iz{a,c,b,b). 
given by 

a = a (3 — b 

and 

a = a /3 = c 



For the step-like initial data (1-19) with a > 
of the Whitham equations (1-18) with g ^ I is 

X 

- = /i3(a,c,-U3,6) , U4 = 6 

Outside this interval, the solution of (1-17) is 

X 

for - < ^i3{a,c,c,b) , 

X 

for - > ^i3{a,c,b,b) . 



Proof. It suffices to show that iJ,3(^a,c,U3,b) is a decreasing function of U3 for 
6 < U3 < c. By (2.10), we have 

dn3{a,c,U3,b) IdXs dq 1 mi 
d^3 =29^i;9^I^ + 2[^^-2(" + " + "^ + ^)]9^- 

The second term is negative because of (2.6) and = 3/8 > 0. The first term is 
also negative: Its first factor is positive in view of (1.13); while its second factor is 

d 1 

= -(a + c+3w3 + &) < 0, 
OU3 4 

for 6 < U3 < c, as we have that a + 6 + 4c < 0. □ 
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Fig. 4.2. Self-Similar solution of the Whitham equations (1.18) with g = 1 and the corre- 
sponding periodic oscillatory solution (1.21) of the mKdV equation with e = 0.018. The initial data 
is given by (1.19) with a = 2, b = —1 and c = of type VI. The oscillations have two distinct 
structures, which are separated by x/t ^ 1.51. 



4.2. Type VI. Here we consider the step initial function (1.19) satisfying > 
a + 56 > -4(c - b) with V{a, c, b) < 0. 

Theorem 4.2. (see Figure ^.2.) For the step-like initial data (1-19) with > 
a + 56 > — 4(c — 6) and V"(a, c, 6) < 0, the solution of the Whitham equations (1.18) 
with g = I is given by 

X X 

ui=a, J = ^i2{a,U2,U3,b) , - = /i3(a, U2, U3, ^) , U4 = b (4.1) 
for ^3{a, —{a + 6)/4, —(a + 6)/4, 6) < x/t < fi3{a, u*** , u*** , 6) and by 

X 

ui = a, M2 = c , y = /i3(a,c,M3,6) , U4 = 6 (4.2) 

for fi^^OjU*** ,u*** ,b) < x/t < ^3{a,c,b,b), where u*** is the unique solution U3 
of 112(0., c,U3,b) — fj,3(a,c,U3,b) in the interval b < U3 < c. Outside the region 
/i3(a, — (a + 6)/4, — (a + 6) /4, 6) < x/t < /is (a, c, 6, 6), the solution of equation (1-17) 
is given by 

X 

a = a, 13 = b for - < ^^3(0, -(a + 6)/4, -(a + 6)/4, 6) , 

and 

X 

a — a , 13 — c for — > ^^{a, c, 6, 6) . 



Proof. We first locate the "leading" edge, i.e., the solution of equation (4.1) at 
U2 — U3. Eliminating x/t from the first two equations of (4.1) yields 

/X2(a, "2, "3, b) - ^i3{a, U2, W3, 6) = . (4.3) 

Since it degenerates at U2 =1*3, we replace (4.3) by 

,x H2{a,U2,U3,b) ~ fi3{a,U2,U3,b) 
G{a,U2,U3,b) := . 77 = . (4.4) 

("2 - ■ti3)v("l ^ '"3)('"2 - Ui)I{a,U2,U3,b) 

In the Appendix, we show that, at the "leading" edge U2 = 1*3, we have 

G{a,U3,us,b)^2{^ + ^)^0 , 
0U2 0U3 
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> 


if 


du2 


du2 




<0 


if '^^ 




du3 



in view of (A. 2), which along with (2.13) gives U2 = U3 = — (a + 6)/4. Having located 
the "leading" edge, we solve equation (4.4) near U2 = W3 = — (a + 6)/4. We use 
formula (A. 3) to obtain 

dG{a,-{a+b)/4,-{a+b)/4,b) ^ dG{a, -{a + 6)/4, -(a + 6)/4, b) ^ ^ 

dU2 dU3 

These show that equation (4.4) gives U2 as a decreasing function of U3 

U2 = B{u3) (4.5) 

in a neighborhood oi U2 = U3 = —{a + b)/4. 

We now extend the solution (4.5) of equation (4.3) as far as possible in the region 
b < U3 < —{a + b)/A<U2 < c. We use formula (2.10) to obtain 

diJ,2 1 d\2 dq 1 

9^ = 2 a^I^ a^I^ + 2 - '('^ + + + 9^11 ' 

P = lP^ + l[Xs-2{a + U2 + us + b)]pl . 
du3 2 du3 du3 2 du3 

In view of (1.13) and (2.6), we have 

fin 

>o , 

<0 . 

OU3 OU3 
We claim that 

|^>0, |^<0 (4.6) 
0U2 0U3 

on the solution of (4.3) in the region & < U3 < —{a + b)/A < U2 < c. To see this, we 
use formula (2.10) to rewrite equation (4.3) as 

l[X2-2{a + U2 + U3 + b)]p- = l[X3-2{a + U2 + U3 + b)]p- . 

L OU2 ^ OU3 

This, together with 

dq dq d^q 1 ) 

du2 du3 du2du3 2 

for U2 > U3, and inequalities (2.6), proves (4.6). 

Hence, the solution (4.5) can be extended as long as 6 < U3 < — (a+6)/4 < U2 < c. 
There are two possibilities; (1) U2 touches c before W3 reaches b and (2) U3 touches b 
before or simultaneously as U2 reaches c. 

Possibility (2) is unattainable. To see this, we use (2.26) to write 

/i2(a, U2, b, b) - /X3(o, U2, b, b) = ^^"^ — V{a, U2,b) , (4.7) 

which is negative for 6 < U2 < c since V{a, U2, b) of (2.27) is an increasing function of 
U2 and since V{a, c, b) < 0. Therefore, U2 will touch c before U3 reaches b. When this 
happens, we have 

/i2(a, c, U3, b) - ii3{a, c, U3, 6) = . (4.8) 
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Lemma 4.3. Equation (4-8) has a simple zero, counting multiplicities, in the 
interval h < < c. Denoting this zero by u*** , then ii2{ci,c,U3,b) — ^^{a, c,U3,b) is 
positive for > u*** and negative for < u*** . 

The proof, which involves formulae (2.24) and (2.25), is rather similar to the proof 
of Lemma 3.19. We will omit it. 

We now continue to prove Theorem 4.2. Having solved equation (4.3) for U2 as 
a decreasing function of U3 for u*** < U3 < —(a + 6)/4, we can then use the middle 
two equations of (4.1) to determine U2 and U3 as functions of x/t in the interval 
/i2(a,-(a + 6)/4, -(a + 5)/4,6) < x/t < ^2(0, c, u***, 6). 

We finally turn to equations (4.2). We want to solve the third equation of (4.2), 
x/t — ii^{a,c,u^,b), for U3 < u*** . It is enough to show that ^,y,{a,c,u^,b) is a 
decreasing function of U3 for M3 < u*** . According to Lemma 4.3, fi2{a,c,U3,b) — 
H^{a,c,u^,b) < for U3 < u*** . Using formula (2.10) for 1^2 and /X3, we have 



d 1 

-[\2~2{a + c + U3 + b)]-^ < -[A3-2(a + c- 

; C/U2 ^ 



U3 



b)] 



dq 
5m.3 



This, together with 



dq 
du2 



dq 
du3 



-(C-U3) > 



for U3 < c, and inequalities (2.6), proves 

dq{a,c, uz,b) 



(9U3 



< 



for M3 < u*** . Hence, 



1 8X3 dq 

2 du^ du3 



+ -[A3 -2(a + c + U3 + 6)] 



d^q 
3 



where we have used inequality (1.13). □ 

4.3. Type VII. Here we consider the step initial function (1.19) satisfying > 
a + 56 > -4(c - b) with V{a, c, b) = 0. 




Fig. 4.3. Self-Similar solution of the Whitham equations (1.18) 5 = 1 and the corresponding 
oscillatory solution (1.21) of the niKdV equation with e = 0.01. The initial data is given by (1.19) 
with a = 1 + b = —I and c = of type VII. 



Theorem 4.4. (see Figure 4-3.) For the step-like initial data (1.19) with > 
a + 5fe > — 4(c — 6) and V{a,c,b) = 0, the solution of the Whitham equations (1.18) 
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with g = 1 is given by 

X X 

ui = a, J = ^i2{a,U2,U3,b) , y = /X3(a, U2, "3, &) • ?i4 = c , 

for fi3{a,—{a + b)/4:,—{a + b)/4:,b) < x/t < iJ,3{a,c,b,b). Outside the region, the 
solution of equations (1-17) is given by 

a = a, P = b, for | < /Z3(a, -(a + 6)/4, -(a + 6)/4, 6) , 



and 



a = a, (3 = c, for - > H3{a,c,b,b) . 



Proof. It suffices to show that U2 and U3 of /Z2(a, M2, U3, 6) — /X3(a, U2, M3, fo) = 
reaches c and b, respectively, simultaneously. To see this, we deduce from equation 
(4.7) that 

H2 (a, c, 6, b) - iiz {a, c, 6, b) = 2(^a + c-2b) ''^ ^'^'^^ 

vanishes when V{a, c, b) = 0. □ 

4.4. Type VIII. Here we consider the step initial function (1.19) satisfying 
> o + 56 > -4(c - b) with V{a, c, b) > 0. 

Theorem 4.5. (see Figure 4-4-) For the step-like initial data (1.19) with > 
a + 56 > — 4(c — 6) and V{a,c,b) > 0, the solution of the Whitham equations (1.18) 
with g =1 is given by 

X X 
Ui=a, - = y^2(a,M2,W3,6) , - = ^3(0, U2, "3, , M4 = 6, 

for /i3(a, — (a+6)/4, — (a+6)/4, 6) < x/t < H3{a,u, u, b), where u is the unique U2-zero 
of the quadratic polynomial V{a, U2, 6) in the interval —(a + 6)/4 < U2 < c. Outside 
the region, the solution of equations (1-17) is divided into the following three regions: 

(1) For x/t < naia, -(a + 6)/4, -{a + 6)/4, 6), 

a = a , (3 = b . 

(2) For ii2{a,u,b,b) <x/t< | (a^ + 2ac+ 50^), 



" = ^=-r+Vl5i-25'^ • 
(3) For x/t>l{a'^ + 2ac+5c^), 

a = a , (3 = c. 

Proof. By the calculation (4.9), when U3 of /Lt2(a, U2, M3, 6) — /Lt3(a, U2, W3, 6) = 
touches 6, the corresponding U2 reaches u, where V{a,ii,b) = 0. Obviously, u < c. 
Hence, equations 

X X 

- = IJ,2{a,U2,U3,b) , - = IJ,3{a,U2,U3,b) 

can be inverted to give U2 and U3 as functions of x/t in the region ij,2{a,—{a + 
b) / A, ~{a + b) / A,b) < x/t < ii2{a,u,u,b). To the right of this region, equations (1.17) 
has a rarefaction wave solution. □ 
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3.2 3.6 4 4.4 3.2 3.6 4 4.4 

x/t X 



Fig. 4.4. Self-Similar solution of the Whitham equations and the corresponding oscillatory 
solution (1.21) of the mKdV equation with e = 0.008. The initial data is given by (1.19) with a = 4, 
b = —1 and c = of type VIII. The solution in the region 4.1 < x/t < 4.2 represents a rarefaction 
wave. 



5. Other Initial Data. We conclude the paper by showing how to handle the 
initial data (1.20). Inequalities (2.14) are replaced by 

d\2 3 Ai — A2 9Ai 

du2 2 Ui — U2 du2 

for U4 < U3 < U2 < wi. Results similar to Lemma 2.2 can then be easily proved. The 
rest of calculations are also similar to those in Sections 3 and 4. 

Appendix A. Leading Edge Calculations. The function /(mi, it2, M3, U4) of 
(2.3) can be written in terms of the complete elliptic integral of the first kind K{s), 
i.e., 

/ - 



^(-Ul - U3)(m2 - -"4) ' 

where 

(Ul - U2)(U3 - M4) 

s 

(Ui - Ui){u2 - Ui) 

Using the derivative formula 

dK{s) _ E{s) - (1 - s)K{s) 
ds ~ 2s(l - s) ' 

where E{s) is the complete elliptic integral of the second kind, we calculate A2 and 
A3 of (2.2) 

Ui — U9 

A2 = 2(71 - 4 



1 - 

A3 = 2cti + 4 



U2—U2 K 
U3 - "4 



^ U2—U4 E_ 

U2 — U2, K 



We then use (2.10) to write /j,2 — A*3 ^ 
-2{u2-U3)K{s) 



ui ~ U2 dq ^ U3 — U4 dq 

{u2 - U3)K - {ui - U3)E du2 {u2 - U3)K - {u2 - Ui)E du3 
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Hence, G(ui, U2, U3, M4) of (4.4) becomes 
ui — U2 dq 



G 



Us - Ui 



dq 



{U2 — U3)K — (Ui - Us)E dU2 (W2 - U3)K — {U2 — U4,)E dus 

We now use the asymptotics of K{s) and E{s) as s is close to 1 

K{s) « 2 ' ^(^) ^l+\{l-s) (log - 1) 



(A.l) 



to calculate the U2 = U3 limit 



G = 2 



dq dq 
du2 duz 



(A.2) 



Finally, we can also use the expression (2.13) and the derivative formula 

dE{s) _ E{s) - K{s) 
ds 2s 

to evaluate the partial derivatives of G in the U2 = W3 limit 

{ui - U4){ui + AU3 + U4) 



dG 

du2 
dG 
duz 



2 + 



= 2- 



4(ui - U3)(U4 - U3) 

(mi — Ui){ui + 4^3 + U4) 

4(mi - W3)(m4 - W3) 



(A.3) 
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